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ABSTRACT 

We present a first analysis of deep 24/im observations with the Spitzer Space Telescope of a sample of nearly 
1500 galaxies in a thin redshift slice, 0.65 < z < 0.75. We combine the infrared data with redshifts, rest-frame 
luminosities, and colors from COMBO-17, and with morphologies from Hubble Space Telescope images collected 
by the GEMS and GOODS projects. To characterize the decline in star-formation rate (SFR) since z ~ 0.7, we 
estimate the total thermal infrared (IR) luminosities, SFRs, and stellar masses for the galaxies in this sample. At 
z ~ 0.7, nearly 40% of intermediate and high-mass galaxies (with stellar masses > 2 x 10 10 M©) are undergoing 
a period of intense star formation above their past-averaged SFR. In contrast, less than 1 % of equally-massive 
galaxies in the local universe have similarly intense star formation activity. Morphologically-undisturbed galaxies 
dominate the total infrared luminosity density and SFR density: at z ~ 0.7, more than half of the intensely star- 
forming galaxies have spiral morphologies, whereas less than ~ 30% are strongly interacting. Thus, a decline in 
major- merger rate is not the underlying cause of the rapid decline in cosmic SFR since z ~ 0.7. Physical properties 
that do not strongly affect galaxy morphology — for example, gas consumption and weak interactions with small 
satellite galaxies — appear to be responsible. 

Subject headings: galaxies: evolution — galaxies: general — galaxies: stellar content — galaxies: fundamental 
parameters — galaxies: interactions — infrared: galaxies 



1. introduction 

It has become clear in the last decade that the average star- 
formation rate (SFR) per unit comoving volume — the so- 
called cosmic SFR — has declined by an order of magnitude 
since z ~ 1 (Lilly et al. 1996; Madau, Pozzetti, & Dickinson 
1998; Hogg et al. 1998; Flores et al. 1999; Haarsma et al. 2000; 
Hopkins 2004). A number of physical processes may contribute 
to this decline: e.g., a declining rate of major galaxy mergers, 
a drop in the rate of minor tidal interactions, or the progres- 
sive consumption of cold gas. Yet, many empirical aspects of 
this declining cosmic SFR are frustratingly unclear. Does the 
star formation in disks evolve slowly, while merger-induced 
starbursts evolve rapidly? Does the decreasing number den- 
sity of galaxies with high SFR reflect a decline in SF in Milky- 
Way mass galaxies or reflect a disappearance of a population of 
bursting dwarf galaxies? In this paper, we present a first basic 
analysis of deep 24/im observations taken by the MIPS Team 
(Rieke et al. 2004; Papovich et al. 2004; Rigby et al. 2004). 
We combine these infrared (IR) data with the COMBO-17 red- 
shift and spectral energy distribution (SED) survey (Wolf et al. 
2004a) and the largest existing Hubble Space Telescope (HST) 
ACS mosaic (GEMS; the Galaxy Evolution from Morphology 
and SEDs; Rix et al. 2004) to study nearly 1500 galaxies with 
0.65 < z < 0.75 in the 30' x 30' extended Chandra Deep Field 
South (CDFS), with the goal of exploring the morphologies, 
stellar masses and SFRs of galaxies, and probing the physics 
driving the decrease in the cosmic SFR to the present day. 



1.1. Understanding the form of the cosmic star formation 

history 

The basic form of the cosmic star-formation history (SFH) 
was initially established using observations of rest-frame ultra- 
violet (UV) luminosity (Lilly et al. 1996; Madau et al. 1996; 
Madau, Pozzetti, & Dickinson 1998; Steidel et al. 1999), which 
in the absence of dust is an excellent indicator of the bolomet- 
ric light output from young, massive stars. As rest-frame UV is 
redshifted into the optical, it becomes observationally straight- 
forward to establish a deep census of unobscured young stars, 
well below the knee of the luminosity function. Yet, studies 
of both local and distant galaxies have shown that only a small 
fraction of the UV photons from ~ L* galaxies escape unhin- 
dered (dwarf galaxies are often optically thin to UV photons, 
but do not significantly contribute to the SFR density; Wang & 
Heckman 1996; Adelberger & Steidel 2000; Bell 2003). Thus, 
the observed drop in UV luminosity density may not track 
quantitatively the drop in cosmic SFR. 

Mindful of the limitations of rest-frame UV, other SFR in- 
dicators have been explored. Rest-frame optical line emission 
reflects the ionizing photon output from very massive young 
stars, in the absence of dust extinction (see Kennicutt 1998, for 
a detailed discussion). Emission line-derived SFRs have been 
studied intensively in the local universe (Gallego et al. 1995; 
Brinchmann et al. 2004), and have been probed out to z ~ 3 
for the very brightest galaxies using rest-frame Oil and Balmer 
line emission (Hogg et al. 1998; Yan et al. 1999; Glazebrook 
et al. 1999; Tresse et al. 2002; Hippelein et al. 2003; Erb et al. 
2003). These studies find a considerable drop in the average 
SFR per unit volume in luminous galaxies since z ~ 1, but ow- 
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ing to the difficulty of measuring Balmer decrements of high 
redshift galaxies, considerable uncertainties remain. 

Radio and X-ray emission from distant galaxies also con- 
strains their SFRs (Condon 1992; Cohen 2003). Observational 
constraints limit the applicability of both indicators to just the 
brightest sources at this stage; a substantial drop in the volume- 
averaged SFR in very luminous sources is again seen since z <~ 1 
(e.g., Haarsma et al. 2000; Cohen 2003). Yet, both X-ray and 
radio luminosities are strongly boosted by the presence of an 
active galactic nucleus (AGN), and the physical origin of both 
is not well-understood. 

One of the most intuitive SFR indicators is the luminosity of 
a galaxy in the thermal IR. For most star-forming galaxies, the 
thermal IR emission is dominated by reprocessed UV photons 
from massive young stars that are dust-enshrouded (see Ken- 
nicutt 1998, for an excellent review). For most galaxies the 
combined IR and UV luminosity is a relatively robust indica- 
tor of the bolometric output from young stars, and therefore a 
good proxy for SFR (e.g., Gordon et al. 2000; Buat et al. 2002; 
Bell 2003), especially if the modest contribution to the IR flux 
from dust-reprocessed light from older stellar populations is ac- 
counted for (typically ~ 30% or less; Misiriotis et al. 2001; 
Bell 2003). Until recently, instrumental sensitivity has been 
a key limitation at these wavelengths; only galaxies with IR 
luminosities > 10 n L Q have been observed at cosmologically- 
interesting redshifts by the Infrared Space Observatory (ISO). 
Through detailed modeling and analysis of these observations, 
it is clear that highly-luminous IR galaxies were much more 
common at z ~ 1 than they are today (e.g., Flores et al. 1999; 
Elbaz et al. 2002; Mann et al. 2002; Pozzi et al. 2004). Observa- 
tions of sub-mm emission from distant, intensely star-forming 
galaxies have extended this picture to z > 2 (e.g., Blain et al. 
1999). Yet, despite this encouraging progress, it is still unclear 
how the evolution of these highly-luminous galaxies is linked 
to the much more numerous fainter galaxies that dominate the 
volume-averaged SFR. 

1.2. Understanding the physical drivers of star formation 

A wide range of processes may contribute to the order-of- 
magnitude drop in cosmic SFR since z ~ 1 : major mergers, tidal 
interactions, gas consumption, ram-pressure stripping, galaxy 
harassment, or the secular evolution of bars. Deciding which 
processes dominate involves a number of observational and in- 
terpretive challenges. As an example, it is challenging to mea- 
sure the major merger rate robustly; close pair statistics suffer 
from contamination, automated asymmetry measures are sen- 
sitive to contamination from irregular galaxies and very close 
projected pairs, and visual morphological classifications are 
subjective, and less sensitive to prograde interactions owing to 
their lack of pronounced tidal tails. 

There is evidence that these different physical processes play 
a role in shaping the SFRs of galaxies. Many merging galax- 
ies form stars intensely (e.g., Sanders & Mirabel 1996). Given 
the observational evidence for a drop in the merger rate since 
z ~ 1 (e.g., Le Fevre et al. 2000; Patton et al. 2002; Conselice 
et al. 2003), it is indeed plausible that the decline in the volume- 
averaged cosmic SFR is driven at least in part by this drop 
in merger rate. Furthermore, minor interactions seem to en- 
hance SFRs, at least in the local universe (e.g., Barton, Geller, 
& Kenyon 2000), plausibly by exciting unstable modes such 
as bars or strong transient spiral arms in the affected galaxies. 
The exhaustion of cold gas could also play a role. Indeed, most 



massive nearby galaxies have low gas fractions < 10% (e.g., 
Kennicutt, Tamblyn, & Congdon 1994; Bell & de Jong 2000), 
and the current SFRs of many local galaxies will exhaust the 
available gas supply in ~ 5 Gyr (e.g., Kennicutt, Tamblyn, & 
Congdon 1994). It is possible that fresh infall or recycled gas 
can help to sustain star formation for a longer time. Yet, it is 
also possible that we live in a special epoch, when the star for- 
mation in massive galaxies is switching off. 

1 .3. The object of this paper 

This paper is one of a series in which we explore different 
aspects of the declining cosmic SFR. Wolf et al. (2005) study 
the contribution of different galaxy morphologies to the rest- 
frame 2800A luminosity density for rest-frame B-band selected 
galaxies in a thin 0.65 < z < 0.75 redshift slice from COMBO- 
17 and GEMS. In this paper, we supplement this optically- 
selected sample with 24/im data from Spitzer, to explore the 
contribution of obscured star formation 1 . Le Floc'h et al. (in 
preparation) explore the IR luminosity function over the range 
< z < 1 using MIPS 24/im data and COMBO-17 redshifts. 
Papovich et al. (in preparation) explore the morphologies of 
the galaxies that dominate the turn-over in the 24/im number 
counts, as a clue to understanding the drivers of IR-obscured 
star formation. 

In this paper, we study the MIPS 24/im properties of a com- 
plete sample of rest-frame B-band selected 0.65 < z < 0.75 
galaxies in the CDFS 2 . Owing to the high sensitivity and rel- 
atively tight point spread function of 24/im Spitzer observa- 
tions, many more star-forming galaxies are detected at 24/im 
than at 70/mi and 160/im. We use IR template spectra, tuned 
to reproduce the whole range of local galaxy IR SEDs, to trans- 
form 24/im luminosity into total IR luminosity. For the smaller 
area imaged by Chandra to date, we incorporate the Chandra 
X-ray data (Alexander et al. 2003) to explore the importance of 
AGN at 24/im (see also Rigby et al. 2004). Our main goals are 
to identify all z ~ 0.7 galaxies with high SFRs, explore their 
morphologies, and characterize the importance of massive star- 
bursting galaxies compared to the present-day universe, as a 
probe of the physics driving the declining cosmic SFR. 

The paper is set out as follows. In §2 we discuss the data. 
In §3 we outline the methods used to estimate total IR flux, 
SFR and stellar mass, and the main sources of uncertainty. In 
§4 we present trends in the IR/UV ratio. In §5 we explore the 
contribution of different morphological types to the total SFR 
at z ~ 0.7. In §6 we discuss intense star formation in massive 
galaxies, and how this changes over the last 7 Gyr. In §7 we 
compare with existing data, explore the contribution of AGN, 
compare with theory, and discuss future improvements in this 
analysis. In §8 we present our conclusions. Appendices A and 
B discuss possible incompleteness in our optically-selected cat- 
alog, and uncertainties in the galaxy morphological classifica- 
tions. Throughout, we express the estimated integrated UV and 
IR luminosities in terms of the bolometric luminosity of the 

'An additional advantage of adopting such a narrow redshift slice for this 
paper is that the effects of differential (1 +z) 4 surface brightness dimming, mor- 
phological ^-corrections, and uncertainties in the conversion of 24/im flux to 
total IR are minimized. 

2 At this stage, we cannot and do not study completely obscured galaxies 
with extremely faint or non-existent UV/optical detections. Preliminary indica- 
tions are that these optically faint IR sources are optically faint largely owing 
to their relatively high redshift (z > 1; e.g., Papovich et al. 2004) rather than 
having a lower redshift and higher obscuration. See Appendix A, where this 
issue is explored in more detail. 
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sun L Q = 3.9 x 10 26 W. We assume n m = 0.3, ft A = 0.7, and 
H = TOkms" 1 Mpc" 1 (Spergel et al. 2003). Estimates of SFR 
and stellar mass assume a universally-applicable Kroupa (2001) 
IMF 3 . 

2. THE DATA 
2.1. MIPS Data 

Spitzer observed a 1° x 0.°5 field around the CDFS in Jan- 
uary and February 2004 as part of the time allocated to the 
Spitzer Guaranteed Time Observers (GTOs). The MIPS 24/im 
data was taken in slow scan-map mode, with individual ex- 
posures of 10 s. We reduced the individual image frames us- 
ing a custom data-analysis tool (DAT) developed by the GTOs 
(Gordon et al. 2005) 4 . The reduced images were corrected for 
geometric distortion and combined to form full mosaics. The 
final mosaic has a pixel scale of 1 "25 pixel" 1 and an image PSF 
FWHM of ~ 6". The Spitzer astrometry is aligned to the ESO 
Imaging Survey (EIS; Arnouts et al. 2001) with a typical ac- 
curacy ^0.5". The overall astrometric solution is then shifted 
by nearly 1" by matching to the much sparser Two Micron All 
Sky Survey catalog (2MASS). This revised astrometric solution 
is matched to < 0.1"with the COMBO-17 coordinate system. 

Source detection and photometry were performed using tech- 
niques described in Papovich et al. (2004); based on the anal- 
ysis in that work, we estimate that our source detection is 80% 
complete at 83 /iJy in the 24/im image. Based on our noise es- 
timates, the photon noise in the 24/mi image is roughly equal 
to the confusion noise at this depth (see, e.g., Dole et al. 2004). 

2.2. Optical Data: COMBO-17 and GEMS 

To date, COMBO-17 has surveyed three disjoint ~ 34' x 33' 
southern and equatorial fields to deep limits in 5 broad and 12 
medium passbands. Using these deep data in conjunction with 
non-evolving galaxy, star, and AGN template spectra, objects 
are classified and redshifts assigned for ~ 99% of the objects 
to a limit of ntR ~ 23.5. Typical galaxy redshift accuracy is 
&/(l + z) ~ 0.02 (Wolf et al. 2004a), allowing construction of 
^0.1 mag accurate rest-frame colors and absolute magnitudes 
(accounting for distance and ^-correction uncertainties). Astro- 
metric accuracy is ~ 0.1". 

To explore galaxy morphology in the rest-frame optical from 
a single observed passband, we study galaxies in one thin red- 
shift slice. Here, we select galaxies from COMBO-17 in the 
CDFS for morphological classification in the narrow interval 
0.65 < z < 0.75 (corresponding to Af ~ 0.5 Gyr, minimizing 
galaxy evolution across this slice). At this redshift, F850LP 
samples roughly rest-frame V-band, and Spitzer's 24/<m pass- 
band samples rest-frame 14/mi, allowing comparison with local 
samples. Furthermore, z ~ 0.7 is close to the median redshift 
of COMBO-17, maximizing the sample size. 

We use F850LP imaging from the GEMS survey (Rix et 
al. 2004) to provide 0.07" resolution rest-frame V-band data 
for 1492 galaxies with 0.65 < z < 0.75. GEMS surveys a 
- 28' x 28' portion of the CDFS in the F606W and F850LP 
passbands to deep limits using the Advanced Camera for Sur- 
veys (Ford et al. 2003) on the HST. The GEMS area is cov- 
ered by a multiple, overlapping image mosaic that includes the 

3 In the case where the stellar IMF varies as a function of time and/or from 
galaxy-to-galaxy, neither SFRs nor stellar masses will be robust except perhaps 
at the order-of-magnitude level. 

4 To date, our preliminary tests indicate that the DAT reduction of MIPS 
24/im data is nearly identical to that provided by the Spitzer Science Center 
pipeline. 



smaller but deeper GOODS area (Giavalisco et al. 2004). One 
orbit per pointing was spent on each passband (63 GEMS tiles 
and 15 GOODS tiles), allowing galaxy detection to a limiting 
surface brightness within the half-light radius of ^f850lp,ab ~ 
24magarcsec" 2 (Rix et al. 2004). At z = 0.7, ACS resolution 
corresponds to ~ 500 pc resolution, roughly equivalent to ~ 1" 
resolution at Coma cluster distances. 

Galaxy classification was carried out by-eye on the z ~ 0.7 
sample by EFB, CW and DHM using the GEMS F850LP imag- 
ing 5 . In what follows, we adopt by-eye classification bins of 
E/S0, Sa-Sd, irregular/compact (Irr), and peculiar/clearly inter- 
acting (Pec/Int). Galaxies were classified on the basis of both 
central light concentration and smoothness. E/S0 galaxies were 
required to have dominant spheroids and smooth light distribu- 
tions. Sa-Sd galaxies were required to have prominent disks 
with signs of ongoing star formation. Galaxies classified as ir- 
regular are chosen to be similar in morphology to low surface 
brightness and Magellanic irregular galaxies in the local uni- 
verse; their irregular light distributions appear to the eye to re- 
sult from stochastic variations and bursts in their star formation 
history rather than from any interaction. In contrast, the pecu- 
liar/clearly interacting designation requires that galaxies have 
signs of tidal features and/or multiple nuclei, with morpholo- 
gies suggestive of a major galaxy merger or merger remnant 6 . 
These classifications are discussed in detail and examples of 
each class shown in Wolf et al. (2005). 

Obviously, these classifications are subjective and depend on 
the depth of the imaging material. In particular, comparison 
with deeper GOODS data for a subset of 290 galaxies in this 
redshift slice has shown that the fraction of faint clearly inter- 
acting galaxies in GEMS is slightly underestimated (Wolf et al. 
2005). While the ideal solution would be to simply consider 
only the GOODS data, number statistics is a key limitation: 
there are ~ 80 galaxies at z ~ 0.7 in the GOODS area that are 
detected at 24/im, compared to the > 400 galaxies detected in 
the whole GEMS area. We therefore choose to study galaxies 
from the whole GEMS area, adjusting the GEMS-only classifi- 
cations as a function of IR luminosity to account for the under- 
estimated faint interaction fraction. We address this topic, and 
the scatter among the three different classifiers where relevant 
in §5 and Appendix B. 

2.3. Cross-correlation of the IR and optical data 

At this stage, we have simply correlated the 24/im imaging 
catalog with the COMBO-17 catalog. The 24/im footprint does 
not exactly match that of COMBO-17; of the 1727 COMBO- 
17 galaxies in our redshift slice, 1436 are covered by the 24/im 
data. A subsample of 1306 of these are covered by the GEMS 
mosaic, and have full morphological, 24/im and redshift in- 
formation. Sources are matched within 2". Roughly 96% of 
24/im sources have only one match in the COMBO-17 catalog. 
The other 4% have two matches in the COMBO-17 catalog; 

5 Automated classification using profile fitting is highly consistent with the 
by-eye classifications that we adopt here, in terms of differentiating between 
bulge-dominated galaxies and disk-dominated galaxies (Bell et al. 2004). Work 
on automated diagnostics of galaxy interactions in GEMS is ongoing, and will 
be discussed in a future work. 

6 In many works galaxies that do not fit onto the Hubble sequence are labeled 
irregular, where no attempt is made to characterize the driver of the irregulari- 
ties. Because we wish to gain insight into the role of interactions in driving the 
star formation activity of galaxies, we attempt, perhaps naively, to distinguish 
between galaxies with clear signs of interaction (Pec/Int) and those without 
(Irr). 
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the closest positional match is taken in these cases 7 . A 24/im 
catalog constructed using optical positions as prior constraints 
would address many of these limitations, but is challenging to 
implement at this stage; we defer a complete analysis of this 
type to a later date. In total, 442 galaxies with 0.65 < z < 0.75 
are detected at 24 fim and have a match in COMBO- 17 (397 
overlap with GEMS); the rest of the COMBO-17 sources with- 
out 24/im detections are assigned 5<r upper limits of 83/dy 8 . 
In appendix A we explore the possibility that our sample selec- 
tion may be biased against highly-obscured IR-luminous galax- 
ies. We find that z ~ 0.7 IR-luminous galaxies should be bright 
enough to be detected and successfully classified by COMBO- 
17, suggesting that the sample selection and main conclusions 
of this paper should be robust to this potential source of incom- 
pleteness. 

3. ANALYSIS 

3.1. Estimating total IR luminosity 

Local IR-luminous galaxies show a tight correlation between 
rest-frame 12-15/im luminosity and total IR luminosity (e.g., 
Spinoglio et al. 1995; Chary & Elbaz 2001 ; Roussel et al. 2001 ; 
Papovich & Bell 2002), with a scatter of - 0.15 dex. Follow- 
ing Papovich & Bell (2002), we use this correlation to con- 
struct total IR luminosity from the observed-frame 24/im data 
(corresponding to rest-frame ~ 14/im). We use the full range 
of Dale et al. (2001) model template spectra to calculate the 
conversion from observed-frame 24/im to 8-1000/im total IR 
luminosity 9 ; we adopt the mean value as our estimate of IR 
luminosity and the RMS as our estimate of conversion uncer- 
tainty. The mean correction factor corresponds approximately 
to Lir ~ 10^„(24/im) to within 20% where Z„(24/im) is the 
monochromatic luminosity at observed-frame 24/im; the tem- 
plate scatter around this mean is larger, ~ 0.3 dex. Different 
choices of template (e.g., Devriendt et al. 1999), or the adop- 
tion of a luminosity-dependent conversion yield similar results. 
The 5cr detection limit of 83/iJy corresponds to limiting IR lu- 
minosity of 6 x 1O 1O L0 or ~ 6M Q yr -1 in terms of an IR-derived 
SFR. 

The conversion from the 24/im photometry to 8-1000/im IR 
luminosity should be accurate a factor of two or better, which is 
adequate for our purposes. There are two types of error. First, 
the conversion depends on fitting templates to the galaxy SEDs. 
There is a total range of a factor of two difference between 
the conversion that would be predicted by various templates 
for galaxies in the appropriate luminosity range of ~ 1O U L0 
(Devriendt et al. 1999; Dale et al. 2001; Chary & Elbaz 2001; 
Lagache, Dole, & Puget 2003; Lagache et al. 2004). Second, 
there is a possibility that the infrared properties of typical high- 
luminosity galaxies differ at z ~ 0.7 from those of the local 
galaxies used in developing the templates. The evidence avail- 
able to us now suggests that the IR SEDs of z > 0.7 seem to be 
adequately spanned by the IR SEDs of local galaxies: the evo- 
lution of the MIR-radio correlation to z ~ 1 (Elbaz et al. 2002; 
Appleton et al. 2004, e.g.,), 15-24/im flux ratios from ISO and 
Spitzer for galaxies at z ~ 1 (Elbaz et al. 2005), and mid-IR 

7 This effect may lead to 1-2% incompleteness if two sources with 24/im 
emission are mistakenly blended into one source 

8 Sources near the map edges where the noise is non-uniform have been 
disregarded in this work. 

'Total 8-1000/im IR luminosities are ~ 0.3 dex higher than the 42.5- 
122.5/tm luminosities defined by Helou et al. (1988), with an obvious dust 
temperature dependence. 



spectra from Spitzer of z ~ 1 ULIRGs (L. Yan, private com- 
munication). This issue will be significantly clarified by future 
work from Spitzer and Herschel. 

3.2. Estimating star formation rates 

We estimate SFRs using the combined UV and IR emission 
of the sample galaxies. This SFR estimator accounts for both 
direct light from young stars from the UV and obscured light 
from the IR, giving a complete census of the bolometric lumi- 
nosity of young stars in the galaxy (e.g., Gordon et al. 2000). 
The primary difficulties affecting this SFR estimator are (1) the 
poorly-constrained fraction of IR light from dust heated by old 
stars, and (2) geometry, in the sense that the IR light is radiated 
isotropically, whereas the UV light escapes along preferred di- 
rections (such as out of the plane of a disk galaxy). Future 
works will be able to address point (1) through more detailed 
modeling of UV/optical/IR SEDs, and will be able to ascertain 
to which extent averaging over inclination angles effectively ad- 
dresses point (2). 

In practice, we estimate the SFR ip using a calibration de- 
rived from the PEGASE stellar population models (see Fioc & 
Rocca-Volmerange 1997, for a description of an earlier version 
of the model), assuming a lOOMyr-old stellar population with 
constant SFR and a Kroupa (2001) IMF: 

7A/(M Q yr- 1 ) = 9.8xlO- 11 x(L IR + 2.2L U v), (1) 
where Li R is the total IR luminosity and Luv = 1-5^^^800 is 
a rough estimate of the total integrated 1216A-3000A UV lu- 
minosity, derived using the 2800A rest-frame luminosity /„.2800 
from COMBO-17. The factor of 1.5 in the 2800A-to-total UV 
conversion accounts for the UV spectral shape of a 100 Myr- 
old population with constant SFR. The SFR is derived assuming 
that Lir reflects the bolometric luminosity of young, completely 
obscured populations, and that Luv reflects the contribution of 
unobscured stars, that must be multiplied by a factor of 2.2 to 
account for light emitted longwards of 3000A and shortwards 
of 1216A by the unobscured young stars. In practice, the SFRs 
are dominated by the IR contribution for the range of SFRs ex- 
plored in this paper. This SFR calibration is derived under iden- 
tical assumptions as the UV and IR SFR calibrations presented 
by Kennicutt (1998). Our SFR calibration is consistent with 
Kennicutt (1998) to better than 30% once different IMFs are 
accounted for. Bell (2003) tested UV+IR-derived SFRs against 
carefully extinction-corrected Ha and radio-derived SFRs for a 
large sample of local star-forming galaxies, finding excellent 
agreement with 0.3 dex scatter and no offset. Accordingly, 
we adopt a systematic error estimate of 0.3 dex (dominated 
by 24/im-to-total conversion error) and 0.4 dex random error 
(with contributions from both 24/im-to-total conversion and 
limitations in UV+IR-derived SFRs). These error estimates are 
supported by comparisons between ISO 15/mi-derived SFRs 
with carefully-constructed Balmer-line SFRs for galaxies at 
< z < 0.8: there was no offset between the SFR scales and 
0.4 dex random scatter for SFRs less than 250M Q yr" 1 (Flores 
et al. 2004). 

3.3. Estimating stellar mass 

To understand how the instantaneous SFR in a galaxy com- 
pares with its mean past rate, it is necessary to estimate the 
existing stellar mass. Under the assumption of a universally- 
applicable stellar IMF there is a tight correlation between rest- 
frame optical color and stellar M/L, that is rather insensitive to 
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FIG. 1. — IR/UV as a function of SFR. Open diamonds show galaxies de- 
tected with MIPS at 24fim, whereas points show 24fim upper limits for galax- 
ies not detected by MIPS; these limits will move towards the lower left with 
a slope of unity as shown by the arrow, and are therefore are consistent with 
the trend observed for the detected galaxies. The dashed lines outline the local 
relation as derived by Adelberger & Steidel (2000). The error bars show our 
uncertainty estimates: 0.3 dex in IR/UV (dominated by 24/im-to-total conver- 
sion uncertainty) and a 0.4 dex SFR uncertainty. 



details of galaxy SFH, dust content and metallicity (e.g., Bell & 
de Jong 2001; Bell et al. 2003; Kauffmann et al. 2003). Based 
on this correlation, we estimate stellar mass using the rest- 
frame B and V band magnitudes: 

log 10 M*/M Q = -0.4(V-4.82) + (-0.628+1.305[fi-V]-0.15), 

(2) 

where the -0.15 dex term converts the stellar masses to a 
Kroupa (2001) IMF 10 . Systematic and random errors in stel- 
lar masses are roughly 0.3 dex, owing primarily to the effects 
of dust, galaxy age and bursts of star formation (Bell & de Jong 
2001; Bell et al. 2003; Franceschini et al. 2003). For stellar 
populations completely dominated by < 1 Gyr-old stars, the lu- 
minosity at a given color is roughly a factor of three higher 
than for galaxies with a wider range of stellar ages. Because 
Eqn. 2 assumes a wide range in stellar age, the stellar masses of 
galaxies with recent very large bursts of star formation will be 
overestimated. 

4. TRENDS IN IR-TO-UV RATIO WITH SFR 

The comparison of UV and 1R luminosities shows that, both 
on a galaxy-to-galaxy basis and globally, the SFR estimates dis- 
cussed below are essentially IR-derived, with a minor < 30% 
contribution from the UV. To assess the relative importance of 
UV and IR light in individual galaxies at z ~ 0.7, Fig. 1 shows 
IR/UV as a function of SFR, as derived from the IR and UV 
luminosities using Eqn. 1 . The range of IR/UV and the trend 

10 More sophisticated stellar mass estimates, based on direct modeling of 
COMBO-17's observed 17-passband SED (Borch et al., in preparation), agree 
with these color-based estimates with no systematic offset and < 40% random 
scatter. 



with total luminosity in the z <~ 0.7 sample is very similar to that 
observed in the local Universe (e.g., Wang & Heckman 1996; 
Adelberger & Steidel 2000; Bell 2003), indicating little or no 
evolution in the relative proportion of obscured star formation 
at least to z ~ 0.7 11 . This is illustrated in Fig. 1, where the 
locus of local star-forming galaxies is also shown 12 . Galaxies 
with SFRs > 6M Q yr" 1 have IR/UV > 10, indicating that un- 
corrected UV SFRs are gross underestimates of the true SFR 
for the most intensely star-forming systems. 

The cumulated total IR luminosity of all detected sources in 
this redshift slice is 7 ± 3 x 10 13 L Q , where the error bar is dom- 
inated by our assumed error in 24/mi to total IR conversion 13 . 
If we assign detections at 83/iJy to all the upper limits, we de- 
rive an upper limit to the total IR luminosity of twice the value 
from the detected galaxies. This result agrees with the number 
counts and our current IR luminosity function estimates, which 
indicate that about 2/3 of the total background at 24/xm is re- 
solved at our detection limit (Papovich et al. 2004, Le Floc'h 
et al., in preparation). We conclude that the total IR luminos- 
ity in this redshift slice is <~ 10 14 Lq. In comparison, the total 
UV luminosity in this redshift slice is 1.5 ±0.3 x 10 13 L Q when 
the luminosity function is extrapolated to zero luminosity (us- 
ing the 2800A luminosity function determined for this redshift 
slice by Wolf et al. 2005) 14 . Thus, at z ~ 0.7, the ratio of IR to 
UV luminosity density is ~ 7^. 

5 . EXPLORING THE RELATIONSHIP BETWEEN STAR FORMATION 
RATE AND MORPHOLOGY 

In this section, we explore which types of galaxies account 
for the bulk of star formation in the z ~ 0.7 universe, when the 
universe was only half of its present age. 

5.1. The relationship between optical properties and star 
formation 

In Fig. 2, we show the rest-frame optical U — V color and V- 
band absolute magnitude of all galaxies in the 0.65 < z < 0.75 
slice, split by visual morphological type. Filled symbols denote 
galaxies detected at 24-pm while open symbols denote galaxies 
for which only 24/im upper limits (83/iJy) were available. In 
both cases, symbol size scales with the logarithm of SFR or its 
upper limit. 

Among visually-classified E/S0 galaxies, those that are blue 
in the optical are also clearly detected at 24/im. Ignoring for the 
moment active nuclei (see §7), we attribute this 24pm flux to 
rapid star formation. In contrast, there are very few galaxies on 

1 1 This is also consistent with a sample of very luminous IR-luminous z ~ 1 
starbursts from a 15/im-selected ISO sample (Aussel et al. 1999, as discussed 
by Adelberger & Steidel). 

12 This has been adapted from their Fig. 11a, where log 10 SFR ~ 
log 10 (Li6oo + A»ldust) _ 9.7, accounting for our adoption of Ho = 
70kms~' Mpc~'. The UV luminosity is estimated by Luv ~ ^/i/ 1600 
(derived using PEGASE for a lOOMyr-old stellar population with constant 
SFR; coincidentally, the UV spectral shape cancels out, leaving 1216A— 3000A 
luminosity equal to monochromatic 1600A luminosity). These estimates 
of SFR and IR/UV should be consistent with our estimates, to within the 
considerable systematic uncertainties. 

13 This is the total luminosity in the volume-limited 0.65 < z < 0.75 sample, 
and so is in essence a luminosity density. We choose to express quantities in 
terms of total luminosities at this time owing to concerns about cosmic variance. 
This topic is considered in more detail by Le Floc'h et al. (in preparation). 

14 It is worth noting that there is a difference in how we define luminosity 
densities, compared with Wolf et al. (2003) and Wolf et al. (2005). These papers 
express luminosity densities in terms of monochromatic solar luminosities at a 
given rest-frame wavelength, whereas we choose to present Luv in terms of 
bolometric solar luminosities. 
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FIG. 2. — The rest-frame U — V color of 1436 optically-selected 0.65 < z < 0.75 galaxies from COMBO-17/GEMS, as a function of visually-classified rest-frame 
V-band galaxy morphology (classifications are adopted at random from three classifiers). Symbol size scales with the logarithm of the estimated SFR: filled circles 
denote galaxies with MIPS 24/im detections (therefore the symbol size denotes SFR estimated using IR+UV) whereas the open circles show upper limits on the 
SFR, derived using the observed UV luminosity and the upper limits on 24/tm flux. 



the red sequence with 24/im detections. Thus, red E/S0 galax- 
ies are IR-faint, which suggests that their red optical colors are 
the result of dominant old stellar populations, not dust-obscured 
star formation. 

Blue visually-classified spiral galaxies (Sa-Sd) are frequently 
detected in the IR. Many optically-red spiral galaxies are IR- 
faint, and are thus red because of their ancient stellar popula- 
tions. Yet, a rather higher fraction of red spirals are detected 
at 24/im compared to red E/S0 galaxies, especially at fainter 
optical luminosities, indicating a more important role for dust- 
reddened spiral galaxies. Morphologically-classified irregular 
galaxies are typically somewhat less luminous, and owing to 
their lower SFRs are less frequently detected by MIPS; how- 
ever, they are qualitatively similar to low-luminosity spirals. 

The vast majority of visually-classified peculiar and clearly 
interacting galaxies are detected at 24/mi; indeed, the few non- 
detections are almost all relatively faint and are not expected 



to be detected at 24/im for plausible UV/IR ratios. Interest- 
ingly, most red Pec/lnt galaxies are strong 24/im sources (in- 
cluding the brightest galaxy at 24/im in the 0.65 < z < 0.75 
slice, which is an interacting galaxy on the blue edge of the 
red sequence). Therefore, we conclude that i) almost all of our 
visually-classified Int/Pec galaxies are intensely star-forming 
at the epoch of observation, and ii) that the optical colors of 
Int/Pec galaxies largely reflect their dust content, rather than 
their SFH. 

In summary, the optical color of morphologically-undisturbed 
galaxies (E/S0 and Sa-Sd galaxies) reflects SFH in a broad 
sense: most red early-type galaxies, and a significant fraction 
of red spiral galaxies, are non-star-forming; whereas the blue 
early and late-type galaxies have important amounts of current 
star formation as probed by the IR. In contrast, clearly inter- 
acting and peculiar galaxies have strong IR-luminous star for- 
mation irrespective of their optical properties (as is expected of 
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Fig. 3 . — The estimated 8-1000/im luminosity function, split by morpholog- 
ical type for 397 galaxies at 0.65 < z < 0.75 galaxies. Only galaxies detected 
at 24/im are shown, and no attempt to extrapolate to lower IR luminosities has 
been made; the sample is grossly incomplete below ~6x 10 10 Lq as denoted 
by the grey dotted line. In each panel, the grey solid histogram shows the 
total IR luminosity function. The shaded area shows the IR luminosity func- 
tion split by galaxy type using GEMS-derived galaxy classifications, where 
the extent of the shaded area shows explicitly the differences in IR luminosity 
function given by the three different classifiers. The black histogram shows 
the IR luminosity function, averaged over the three different classifiers and 
corrected to reproduce the increased fraction of clearly interacting galaxies 
seen in GOODS-depth data, as described in Appendix B. 



very dusty star-forming systems). 

5.2. The contribution of different galaxy types to the cosmic 

SFR 

An important question is the relative contribution of galax- 
ies with different morphological types to the total SFR density 
at z ~ 0.7. The cosmic-averaged SFR has dropped by a factor 
of three since z ~ 0.7 as measured in the IR (e.g., Flores et al. 
1999; Pozzi et al. 2004), radio (e.g., Haarsma et al. 2000), Ha 
(e.g., Tresse et al. 2002; Perez-Gonzalez et al. 2003), or UV 
(e.g., Lilly et al. 1996; Schiminovich et al. 2005). Thus, an un- 
derstanding of the relative contribution of different galaxy types 
to the cosmic SFR at z ~ 0.7 will constrain the types of galaxy 
that need to 'switch off their star formation activity between 
z ~ 0.7 and the present day, giving insight into the physical 
mechanisms driving this evolution. 

In Fig. 3, we show the distribution of estimated IR luminos- 
ity as a function of galaxy type for 24/im-detected galaxies at 
z ~ 0.7; the top axes show the equivalent SFRs. We choose not 
to normalize the luminosity function to the observed volume 
at this stage owing to concerns about cosmic variance (e.g., 
Somerville et al. 2004); a thorough analysis of evolution in 
the IR luminosity function is presented in a companion paper 
by Le Floc'h et al. (in preparation). We do not extrapolate 
the IR luminosity function below the 24/im detection limit at 
~6xl0 10 L Q . 



The IR luminosity density for 10.7 < log 10 LiR/L Q < 11.5 
is dominated by spiral galaxies, with an important contribution 
from clearly interacting galaxies. Morphologically-classified 
E/S0 galaxies and irregulars contribute, but play lesser roles. 
Clearly interacting galaxies make an increasingly important 
contribution at log 10 Lr R /L Q > 11.3. This behavior resem- 
bles that in the local universe, where almost all galaxies with 
log lo Li R /L > 11.8 show clear signs of interaction (Sanders 
& Mirabel 1996). Integrating over the whole distribution of 
galaxies detected at 24/im, we find that the fractional contribu- 
tions of spiral, E/S0, Pec/Int and irregular galaxies to the total 
integrated IR luminosity density are 0.51 : 0.10 : 0.28 : 0.11. 

To probe the uncertainties in the relative contributions to the 
SFR density from different galaxy types, we explore two lim- 
iting cases: (1) all 24/im non-detections have no IR flux (i.e., 
their SFR is measured purely from the UV); or (2) all 24/im 
non-detections emit at the detection upper limit at this wave- 
length of 83/iJy. In the first case, the split of SFR is (Spiral : 
E/S0 : Pec/Int : Irr) 0.51 : 0.10 : 0.26 : 0.13, very similar to the 
split in terms of observed IR flux (because the UV is a relatively 
small contribution to the total SFR density). In the second case, 
where limits are all taken to be marginal detections, the split 
is 0.47 : 0.15 : 0.18 : 0.20. Because of the preponderance of 
E/S0 and Irr non-detections, the relative importance of these 
classes is increased (arguably, rather artificially). In neither of 
these cases are the contributions from optically-faint galaxies 
(i.e., those galaxies fainter than M v > -19) accounted for. Be- 
cause the integral over a luminosity function is dominated by 
~ L* galaxies, our limitation to relatively bright sources does 
not prevent us from qualitatively constraining the contribution 
of different galaxy types. Nonetheless, it is clear that the de- 
tailed form of the distribution of SFRs to lower limits will im- 
portantly affect the quantitative contributions of the different 
galaxy types, and is arguably the largest source of systematic 
uncertainty in this paper. 

To summarize, the SFR density at z ~ 0.7 is dominated 
by morphologically-normal spiral, E/S0, and irregular galaxies 
(> 70%). Clearly interacting galaxies account for < 30%. Our 
result is based on a combination of UV and IR luminosities, 
but is dominated by IR-luminous objects. Wolf et al. (2005) 
have examined the UV-luminous objects alone and also find 
that most of the rest-frame UV flux density at z ~ 0.7 is in 
morphologically-undisturbed galaxies. 

These results have important implications for the physical 
mechanisms which drive the factor-of-three decline in cosmic- 
averaged comoving SFR density from z ~ 0.7 to the present 
day. Because clearly interacting galaxies — with morphologies 
suggestive of major galaxy mergers or merger remnants 15 — 
contribute much less than half of the z ~ 0.7 SFR density, the 
drastic decline in cosmic SFR is only partially driven by the 
decreasing frequency of major galaxy mergers at the present 
day (e.g., Le Fevre et al. 2000; Patton et al. 2002). Instead, the 
dominant driver is a strong decrease in SFR in morphologically- 
undisturbed spiral galaxies. It is possible that the decreasing 
star-formation activity of spiral galaxies is driven primarily by 
gas supply and consumption in a purely quiescent mode of star 
formation. However, it is also possible that even in spiral galax- 
ies much star-formation activity is triggered by relatively minor 

15 It is worth noting that relatively advanced merger remnants are also classi- 
fied as clear interactions, until their tidal features have dissipated. Owing to the 
long dynamical times at large galactic radii, we expect that major mergers are 
classified as clearly interacting well after the major final burst of star formation 
has ceased. 
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tidal interactions, that function to trigger internal instabilities 
such as spiral arms or bars, enhancing their SFR. Our observa- 
tions cannot sensitively test these competing hypotheses at this 
stage: deeper redshift surveys and exploration of their gas con- 
tent through, e.g., CO mass determination from the Atacama 
Large Millimeter Array will start to allow us to test the relative 
importance of these physical processes. 

6. EXPLORING THE DEMISE OF STAR FORMATION IN MASSIVE 
GALAXIES 

We turn now to explore the relationship between SFR and 
stellar mass at z ~ 0.7, comparing it to z ~ 0. We choose to ex- 
plore the 'specific SFR' ; i.e., the SFR per unit stellar mass. This 
is closely related to the more model-dependent birthrate param- 
eter, b = SFR/ (SFR), where (SFR) is the past-averaged SFR. 
The past-averaged SFR is estimated by dividing the stellar mass 
by the elapsed time since star formation began in a given galaxy 
(see, e.g., Kennicutt, Tamblyn, & Congdon 1994), adjusting for 
the mass lost by the stellar population as it ages. Roughly 50% 
of initially-formed stellar mass is lost within one billion years 
for a Kroupa (2001) IMF; therefore b ~ SFR/(2M»/f. s/ ). In 
estimating t s f, we simplistically assume the epoch of first star 
formation is at Zf = 4 for all galaxies; using a present-day age 
of the Universe of 13.5 Gyr, this corresponds to t s f ~ 12Gyr 
at z ~ and t s j ~ 6 Gyr at z ~ 0.7. Despite the extra model 
dependence introduced by having to assume a timescale when 
deriving b, birthrate has the attractive feature that it compares 
the observed SFR with the SFR required to build up the existing 
stellar mass within a Hubble time (at that epoch). 

In panel a) of Fig. 4, we show specific SFR as a function of 
stellar mass for the 1436 0.65 < z < 0.75 galaxies in our sam- 
ple. Owing to the detection limit of the 24^im observations, it is 
impossible to detect galaxies forming stars at anything less than 
a constant rate (SFR/M, < 0.4, or equivalently b < 1) for stellar 
masses less than ~2x 10 10 M Q . Yet, despite this rather strin- 
gent limitation, large numbers of galaxies have been detected at 
24/im, with typical specific SFRs of ^1/Gyr (i.e., if the current 
SFR were to continue at its present rate the stellar mass of the 
galaxy would double within 2 Gyr, remembering that 50% of 
the mass initially in stars is returned to the gas phase by stellar 
winds and supernovae). 

Given the existence of an important population of intensely 
star-forming galaxies at z ~ 0.7, it is instructive to explore 
how this compares with the properties of star-forming galax- 
ies in the local universe. We have therefore constructed a lo- 
cal 'control' sample of galaxies, which we show in panel b) 
of Fig. 4. We use the NASA/IPAC Extragalactic Database 
(NED) to select a pseudo volume-limited sample in the range 
1500 < cz/kms" 1 < 3000 with 2MASS K-band magnitude cut 
K < 12, and galactic latitude b > 30deg. This sample is com- 
plete for galaxies above Mk < -19.6. Where they were defined, 
IR fluxes were taken also from NED; total 8-1000/im IR fluxes 
were estimated following Appendix A of Bell (2003), and are 
accurate to <~ 30%. Upper limits on IR flux were roughly esti- 
mated assuming total IR fluxes just at the IRAS detection limit 
of ~ 3 x 10 _11 ergscm _2 s _1 (the conclusions do not depend on 
the choice of this limiting flux). A legitimate and unavoid- 
able concern is the inhomogeneous completeness properties of 
NED redshifts. For the purposes of volume estimation, we es- 
timate an overall NED redshift completeness from the number 
of K < 12, b > 30deg galaxies with redshifts divided by the 
number of K < 12, b > 30deg galaxies both with and without 
redshifts: 14350/23658 (or 61%). Thus, the estimated total vol- 
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FIG. 5. — The cumulative distribution of the birthrates of massive galaxies 
with M* > 2 X 1O 1O M0. The grey shaded region denotes the area over which 
upper limits on birthrates contribute to the distribution (i.e., the shape of the 
cumulative distribution in the grey shaded area is strongly biased towards high 
birthrates). 



ume probed is 71500 x 0.61Mpc 3 . The total sample is 2177 
galaxies, 1089 with IR detections and 1088 with only IR up- 
per limits. In fact, because of large-scale overdensities in this 
volume, there are <~ 3 times more galaxies in this volume than 
one would expect on the basis of, e.g., the stellar mass func- 
tions or TT-band luminosity functions of Cole et al. (2001) or 
Bell et al. (2003), taking into account all the relevant selection 
limits. To reflect this overdensity, we assign an 'effective vol- 
ume' to the sample of ~ 1.2 x 10 5 Mpc 3 ; this 'effective volume' 
matches the actual volume in the 0.65 < z < 0.75 slice. Accord- 
ingly, we use this control sample unaltered, with the rationale 
that in the two panels of Fig. 4 we are comparing galaxies from 
the same effective volumes at both cosmic epochs. None of our 
conclusions are affected by the detailed choice of volume or, in- 
deed, comparison sample. Different NED-derived samples can 
be chosen, or a sample of GEMS/MIPS galaxies with redshifts 
0.1 < z < 0.2 can be used as a control, with unchanged results 
but a reduced number of comparison galaxies. 

Stellar masses and SFRs for the local comparison sample 
were estimated assuming a Kroupa IMF. Stellar masses were 
estimated directly from the /f-band absolute magnitudes (where 
Hubble-flow distances are used, for simplicity) assuming a sin- 
gle TT-band stellar M/L of 0.6 in solar units, appropriate for a 
Kroupa IMF (following Bell et al. 2003; Cole et al. 2001). SFRs 
were estimated from total IR flux following equation 4 of Bell 
(2003) adjusted to a Kroupa IMF, and are not corrected for any 
(typically modest; < 30%) contribution to the dust heating from 
old stellar populations (like the SFR estimates constructed us- 
ing the 24^m data). The stellar masses and SFR estimates are 
accurate to <~ 0.3 dex, assuming a universally-applicable stellar 
IMF, and are directly comparable with stellar masses and SFRs 
for the z ~ 0.7 redshift slice. 
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FIG. 4. — Specific SFR as a function of galaxy stellar mass, in solar units, a) SFR per Gyr per unit stellar mass as a function of stellar mass for the z ~ 0.7 galaxy 
sample. Plotted on the right-hand axis is the birthrate parameter b = SFR/(SFR). Galaxies detected at 24^m are shown with open circles; non-detections have 
upper limits shown by points, b) SFR per Gyr per unit stellar mass as a function of stellar mass for the local comparison sample. Again, the birthrate parameter is 
plotted on the right-hand axis for comparison (the offset between panels a and b is because of the difference in the age of the universe at the two redshift intervals of 
interest). IR detections are shown as open circles; non-detections are denoted by points. The dotted line shows the minimum mass for which the sample is complete. 
In both panels, the solid grey line shows a birthrate of unity; i.e., a constant SFR through all of cosmic history. The grey box shows the position of intermediate- and 
high-mass galaxies with ongoing bursts of significant star formation (b > 1). Both datasets sample ~ 1.2 X 10 5 Mpc 3 , therefore direct comparison of both plots is 
fair, remembering the influence of the detection limits. 



We compare the distribution of specific SFRs and stellar 
masses in the z ~ 0.7 sample and the local control sample in 
Fig. 4; the cumulative histogram of birthrates is given in Fig. 5. 
Both datasets effectively sample <~ 1 .2 x 10 5 Mpc 3 , therefore in 
the regions where both datasets are complete it is fair to com- 
pare the relative distributions of datapoints. The local sample 
is characterized by a broad distribution of star-forming galax- 
ies with specific SFRs between 0.03 and 0.2 (0.1 < b < 1; i.e., 
slowly declining SFR through cosmic history). Furthermore, 
where the detection limits permit, one observes a sizeable pop- 
ulation of galaxies with very low specific SFRs < 0.01, indica- 
tive of a largely non-star-forming subpopulation of (primarily 
early-type) galaxies. There are very few galaxies with bursts of 
star formation that are adding significantly to their existing stel- 
lar mass (i.e., birthrates b > 1). In particular, at stellar masses 
in excess of 2 x 10 10 M Q , where the sample is essentially com- 
plete, 5/699 galaxies (< 1%) of galaxies have b > 1 (the grey 
shaded region) 16 . While there are significant systematic un- 
certainties in both SFR and stellar mass, it is impossible to ac- 
commodate a large population of starbursting intermediate- and 
high-mass galaxies in the local universe. This result is insensi- 
tive to the adoption of a variety of plausible prescriptions for 
SFR and stellar masses, and indeed the rarity of galaxies with 
b > 1 at M* > 2 x 1O 1O M is apparent from Fig. 24 of Brinch- 
mann et al. (2004), who analyse accurate emission line-derived 
SFRs from the Sloan Digital Sky Survey. 

In stark contrast, turning to the z ~ 0.7 sample, there are 
a large number of massive galaxies with intense star forma- 
tion. Quantitatively, at stellar masses in excess of 2 x 1O 1O M , 
120/3 1 1 (39%) of galaxies are forming stars at a rate faster than 
their past-averaged rate (b > 1; the grey shaded region in Fig. 
4) 17 . While there are always significant systematic uncertain- 

16 Galaxies with M* > 2 x 1O 1O M0 contain roughly 70% of present-day stel- 
lar mass (Bell et al. 2003). 

17 If starbursts were defined by high specific SFRs instead, the difference 



ties, including estimation of total IR luminosity, SFRs and stel- 
lar masses, these differences are large enough to be robust to 
systematic uncertainties of greater than a factor of three 18 . In 
addition, where possible we have minimized the opportunities 
for systematic error by using IR-derived SFRs in both cases, 
and equivalent SFR and stellar mass calibrations. Indeed, adop- 
tion of the 0. 1 < z < 0.2 GEMS/MIPS sample as the 'local' con- 
trol sample yields only 1 galaxy with b > 1 out of 41 galaxies 
with M, >2x 1O 1O M . In this case, identical methodologies 
have been used for the 0. 1 < z < 0.2 and 0.65 < z < 0.75 sam- 
ples. 

Recalling the conclusions reached in §5, it is worth briefly 
commenting on the visual morphologies of galaxies with M* > 
2 x 1O 1O M . 70/124 massive spiral (Sa-Sd) galaxies, 4/5 mas- 
sive irregular galaxies, and 18/23 Pec/Int galaxies have b > 1, 
i.e. are undergoing significant bursts of star formation. In con- 
trast, only 20/138 massive E/S0 galaxies are forming stars in- 
tensely 19 . Therefore, roughly 1/2 of spiral, irregular and clearly 
interacting galaxies are undergoing starbursts at z ~ 0.7. Fur- 
thermore, the starbursting population is dominated by spiral 
galaxies rather than Pec/Int galaxies. This reinforces the con- 
clusions of §5; even amongst the most intensely star-forming 
z ~ 0.7 galaxies, the vast majority are not forming stars in- 
tensely because of a recent major merger. 

Taken together, it is fair to conclude that while a large frac- 
tion of intermediate- and high-mass galaxies are forming stars 

between z ~ 0.7 and the present day would be even more striking. 

18 It is likely that our neglect of heating of dust by old stellar populations 
affects galaxies with low b values preferentially, artificially inflating their esti- 
mates of SFR and therefore b. From inspection of Fig. 4 this is likely to affect 
the most massive galaxies the most severely, and will affect the low redshift 
sample more acutely than the z ~ 0.7 sample. Thus, if the heating of dust by 
old stars were accounted for by detailed SED modeling, it is likely that the 
conclusions would remain unchanged, or further strengthen. 

I9 The total is 290 galaxies, which is lower than the COMBO-17+MIPS num- 
ber of 31 1 galaxies owing to the smaller areal coverage of GEMS. 
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intensely 7 Gyr ago, almost none are at the present epoch. Fur- 
thermore, these intensely star-forming galaxies are primarily 
morphologically undisturbed, indicating that major mergers are 
not responsible for the bulk of this intense star formation. In- 
stead, physical properties that do not strongly affect galaxy 
morphology — for example, a larger abundance of gas at ear- 
lier times or weak interactions with small satellite galaxies — 
appear to be responsible for the high SFRs. 

7. DISCUSSION 
7.1. Comparison with previous results 

Our results confirm and extend the previously-observed rapid 
evolution in the number of IR-luminous galaxies at intermedi- 
ate and high redshift (e.g., Elbaz et al. 1999; Chary & Elbaz 
2001; Papovich et al. 2004), and with the drastic decline in the 
space density of intensely star-forming galaxies as probed with 
a variety of different SFR indicators (e.g., Cowie et al. 1996; 
Wolf et al. 2003; Bauer et al. 2005). Furthermore, we have 
been able to estimate the stellar masses of these intensely star- 
forming systems, finding that almost all of the intensely star- 
forming galaxies with L 1R > lO n L (luminous infrared galax- 
ies; LIRGs) seen at intermediate redshift have stellar masses 
in excess of ~ 1O 1O M (Fig. 4). This is in agreement with, 
e.g., Franceschini et al. (2003) and Zheng et al. (2004), who 
used optical and NIR data to estimate the stellar masses of ISO- 
CAM 15yum-selected LIRGs, finding stellar masses in excess of 
> 1O 1O M , and Le Floc'h et al. (2004), who explore the proper- 
ties of 26 1 < z < 2.5 IR-luminous sources with Spitzer, finding 
significant stellar masses for their sample of LIRGs. 

Our conclusion that a decline in major interaction rate is not 
the primary cause of the decline in cosmic SFR since z ~ 0.7 
also confirms and extends previous observations. Follow-up 
observations of fields targeted by ISO tentatively suggested that 
the rapidly-evolving IR-luminous population may not be com- 
posed exclusively of strongly-interacting galaxies. Flores et al. 
(1999) found that only 6/16 galaxies in a sample with ISO and 
HST data had strongly disturbed morphologies, while the rest 
were disk-dominated, E/S0, or unresolved. Zheng et al. (2004) 
extended their work to a sample of 36 ISO-selected LIRGs at 
z < 1.2, finding that 17% of their sample were obvious merg- 
ers, and placing an upper limit of 58% on the merger fraction. 
Our study places these indications on firm ground, pushing to 
fainter limits of ~ 6 x 1O 1O L with an order of magnitude larger 
sample of galaxies. 

7.2. How much dust can be heated by AGN at z ~ 0.7 ? 

It is worth discussing the contribution of active galactic nu- 
clei (AGN) to the 24pm total luminosity density. To aid in iden- 
tifying AGN, we use the Chandra data for the CDFS to select 
AGN candidates. The Chandra footprint on the CDFS (Alexan- 
der et al. 2003) is considerably smaller than the GEMS and 
MIPS coverage of the field, therefore in this section we consider 
only galaxies falling completely within the Chandra coverage. 
This area contains 3122 galaxies with COMBO- 17 photomet- 
ric redshifts; the small subsample classified by COMBO-17 as 
galaxies with 0.65 < z < 0.75 and detected by Chandra com- 
prises 22 galaxies; 0.5-8 keV source fluxes and errors are taken 
from Alexander et al. (2003). 

While highly-luminous X-ray sources are invariably AGN, 
at fainter limits a contribution from star-forming galaxies is ex- 
pected. To test this possibility, we compare 24^m-derived IR 
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Fig. 6. — A comparison of 0.5-8 keV observed frame Chandra X-ray lumi- 
nosities with 24/im-derived total IR luminosities for the 22 0.65 < z < 0.75 
Chandra- selected galaxies with optical detections. Filled circles denote MIPS 
detections, whereas open circles with limits denote MIPS upper limits only. 
The IR luminosity errors include uncertainty in conversion to total IR. The 
shaded area shows the expected locus of star-forming galaxies. 

luminosities and X-ray luminosities for these 22 z ~ 0.7 galax- 
ies in Fig. 6. The shaded area shows the expected locus of star- 
forming galaxies, derived from the empirical relationship be- 
tween X-ray and fe-corrected 1 .4GHz radio luminosities (Fig. 1 
of Cohen 2003) coupled with the radio SFR calibration of Bell 
et al. (2003), accounting for the 0.3 dex empirical scatter in the 
correlation 20 . It is clear that only one galaxy has X-ray emis- 
sion that is consistent with a star formation origin; the others 
have excess X-ray emission indicative of an AGN contribution. 

Considering the 21 AGN-dominated X-ray sources, we find 
a total IR luminosity between 3.4 x 1O 12 L and 4 x 1O 12 L , 
where in the first instance we ignore upper limits and in the 
second instance we take upper limits to represent marginal de- 
tections. The total IR luminosity in all detected sources in the 
Chandra footprint is 2.6 x 1O 13 L , increasing to 5.4 x 1O 13 L if 
all non-detections are taken as marginal detections at the 83/iJy 
limit. Accordingly, we estimate that < 15% of the total IR 
0.65 < z < 0.75 luminosity density is in sources with significant 
AGN emission. Clearly, the contribution to the IR luminosity 
density of AGN-heated dust may be substantially lower, as the 
IR luminosity in these galaxies may come from star formation. 
On the other hand, we have not accounted for any contribution 
of AGN to the IR luminosities of galaxies not detected by Chan- 
dra; bearing in mind that the Chandra observations reach deep 
enough to start probing star-formation-powered X-rays from 
galaxies with Ltr ~ 1O U L , one hopes that the AGN in non- 
AGN dominated galaxies will be weak and will contribute rel- 
atively little to the total IR. This conclusion is consistent with 

20 This locus does not account for the mild X-ray ^-corrections in the ob- 
served frame 0.5-8 keV X-ray luminosity, but was defined for galaxies with 
0.4 < z < 1 .3 with a median redshift of 0.76 and so will be accurate to ~ 0. 1 dex 
(Cohen 2003). It is worth noting that the use of the SFR calibration of Cohen 
(2003) results in a similar locus to within 0.1 dex. 
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recent ISO and Spitzer results, finding 10-20% contributions at 
15/im and 24 /im respectively (Fadda et al. 2002; Franceschini 
et al. 2005), and with Silva, Maiolino, & Granato (2004) who 
model the AGN contribution to the IR background, suggesting 
that 10-20% of the IR background at mid-IR wavelengths is 
from AGN and their host galaxies (the AGN alone account for 
5% of the IR background in this model). 

7.3. Comparison with Theoretical Expectations 

Gas-dynamical processes, such as star formation or feed- 
back, are currently very challenging to include in models of 
galaxy formation and evolution owing largely to important res- 
olution limitations. As a result of the complexity of these pro- 
cesses, it has been impossible to converge on a robust descrip- 
tion of star formation; instead, a wide range of empirically- 
motivated prescriptions on kpc or larger scales have been adopted 
(e.g., Somerville & Primack 1999; Springel & Hernquist 2003a; 
Barnes 2004). Current models tend to reproduce the cosmic 
SFH reasonably well (e.g., Cole et al. 2000; Somerville et al. 
2001; Springel & Hernquist 2003b). A few models have started 
to incorporate detailed dust prescriptions, and can model the 
evolution of IR-bright galaxies (e.g., Granato et al. 2000; Bal- 
land et al. 2003). While the progress is very encouraging, de- 
tailed observations such as the 24/im number counts are rather 
poorly reproduced by the models (Papovich et al. 2004), re- 
flecting the difficulty of the 'gastrophysics' that must be tack- 
led to properly reproduce the full richness of the observational 
datasets. For this reason, we focus in this section on qualitative 
features of the models, rather than a detailed comparison of IR 
luminosity functions, or the distribution of stellar masses and 
morphologies as a function of IR luminosity. 

Somerville et al. (2001) present a 'collisional starburst' model, 
in which major (with mass ratios <4:1) and minor (with mass 
ratios between 4: 1 and 10: 1) mergers trigger an episode of rapid 
star formation, with a star formation efficiency that increases to- 
wards lower mass ratio (i.e., is higher for major mergers). At 
z ~ 0, Somerville et al. (2001) predict that 35% of star forma- 
tion is triggered by minor mergers, and 2.5% of star formation 
is triggered by major mergers. At z ~ 0.7, 50% of star forma- 
tion is triggered by minor mergers and 3.5% by major merg- 
ers. Thus, they predicted that the steep drop in the cosmic SFR 
since z ~ 1 is not caused by changes in the major merger rate; 
however, in their model a significant part of the drop reflects a 
declining rate of minor mergers. There are a number of impor- 
tant uncertainties: a lowered star formation efficiency in minor 
interactions would shift more star formation into major merg- 
ers because more gas would be available to form stars in major 
mergers, tidal interactions (e.g., flybys) are not included in the 
model, untriggered star formation in a galaxy with triggered 
star formation is not included in the census of triggered star 
formation, and our observational census of 'clearly interacting' 
galaxies may well include a significant fraction of minor merg- 
ers and tidal interactions. 

Nonetheless, it is encouraging that the observations and at 
least this model are reasonably consistent in the sense that ma- 
jor galaxy mergers do not drive the declining cosmic SFR from 
z ~ 1 to the present day. At present, models cannot clearly pre- 
dict which processes dominate: while minor mergers and tidal 
interactions may dominate (Somerville et al. 2001), it is also 
possible that a dwindling gas supply could dominate (Cole et 
al. 2000; Somerville et al. 2001), depending on the efficiency 
of triggered star formation in the real universe. 



7.4. Improvements 

There are a number of areas that need to be improved to more 
fully understand these results, and to illuminate the physics 
driving the cosmic SFH. We will mention a few specific ex- 
amples here. A more complete, ground-based and IRAC near- 
IR-selected photometric and/or spectroscopic sample covering 
comparable or larger areas will be required to push our under- 
standing past z ~ 1, into the regime where we observe rapid 
cosmic star formation. A wide range of morphological distur- 
bance diagnostics must be developed and robustly tested (see, 
e.g., Conselice et al. 2003; Lotz, Primack, & Madau 2004, for 
some promising examples of this kind of technique). Total Lir 
and SFR calibrations must be tested thoroughly using IR spec- 
tral data coupled with 70^m and 160/mi imaging data where 
available. This must be accompanied by a fuller analysis of lo- 
cal, well-studied control samples to explore the fraction of Lir 
coming from dust heated by optical light from old stellar popu- 
lations, and possible metallicity dependences in IR SEDs. The 
IR properties of high redshift AGN must be better constrained, 
including constraints on the fraction of light emitted by AGN- 
heated dust compared with dust heated by associated star for- 
mation in the AGN host. Yet, the two main conclusions of this 
paper — the rapid demise of intermediate- and high-mass star- 
forming galaxies between z ~ 0.7 and the present-day, and the 
limited importance of major mergers in driving the declining 
cosmic SFR — are likely to hold even after these studies are 
completed. 

8. CONCLUSIONS 

In this paper, we have explored the 24/mi-derived IR prop- 
erties of an optically-selected sample of 1436 galaxies with 
0.65 < z < 0.75 in the CDFS. Optical IDs, photometric red- 
shifts and morphologies were obtained from the COMBO- 17 
and GEMS surveys. From examination of local LIRGs and 
ULIRGs, and exploration of COMBO-17's photometric red- 
shift accuracy as a function of IR luminosity, we argue that the 
optical selection should not introduce any special bias against 
intensely star-forming obscured galaxies. Under the assump- 
tion that the mid-IR properties of z ~ 0.7 galaxies are spanned 
by calibrating samples in the local universe, we have derived 
total IR luminosities for the 442 galaxies detected at better than 
5a at 24/im. We used these data, in conjunction with the optical 
data, to estimate SFRs and stellar masses. 

We have found the following: 

• Galaxies with X-ray-luminous AGN contribute < 15% 
towards the integrated z ~ 0.7 IR luminosity density. 

• The IR-to-UV ratio is a strong, increasing, function 
of total SFR; this relation shows a > 0.5 dex scatter. 
There is no evidence for evolution of this relationship 
over the last 7 Gyr. 

• Morphologically-normal galaxies show a strong 
relationship between their optical and IR properties. 
Optically blue early-type and spiral galaxies are 
IR-bright. Optically red early-type galaxies and a 
significant fraction of red spiral galaxies are not 
detected at 24/tm Most galaxies on the red sequence 
are genuinely old, red and dead. 

• Dust strongly shapes the optical properties of 
morphologically-peculiar galaxies. The majority of 
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obviously interacting galaxies are detected at the star 
formation-sensitive 24-fim band, despite showing the 
full range of optical colors from very blue through to 
very red. 

• Clearly interacting galaxies with morphologies 
suggestive of major galaxy mergers contribute at most 
30% of the integrated IR luminosity density at z ~ 0.7. 
Bearing in mind that there is a factor of three reduction 
in cosmic IR luminosity density between z ~ 0.7 and 
the present day, this suggests that a declining major 
merger rate can only form a small part of this drop in 
cosmic SFR. 

• Morphologically-normal galaxies (spirals, ellipticals 
and Magellanic irregulars) form more than 70% of 
the integrated IR luminosity density in the z ~ 0.7 
slice. This implies that physical processes which do 
not strongly affect galaxy morphology, such as gas 
consumption and weak tidal interactions with small 
satellite galaxies, are likely responsible for the drastic 
decline in cosmic SFR since z ~ 0.7. 

• At z ~ 0.7, nearly 40% of intermediate- and high-mass 
galaxies are actively starbursting (i.e., they have 

M, >2x 1O 1O M and b > 1). In contrast, only 1% of 
z ~ galaxies with > 2 x 1O 1O M are undergoing 
such active star formation. The declining SFR since 
z ~ 0.7 is driven by drastic evolution in the SFRs 
of intermediate- and high-mass galaxies, rather than 
starbursting dwarf galaxies. 
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FIG. A7. — Rest-frame U — V colors as a function of absolute magnitude in 
V-band for a sample of ULIRGs and LIRGs in the local universe (open dia- 
monds). Also plotted are visually-classified peculiar and interacting galaxies 
from the z ~ 0.7 slice, reproduced from Fig. 2. 



A. ARE MOST STARBURSTS IN THE COMBO- 17 SAMPLE? 

A possible source of bias could be that the most heavily- 
obscured z ~ 0.7 star-forming galaxies may be too optically 
faint to be reliably picked up and/or correctly classified by 
COMBO-17. Indeed, there are many MIPS sources without 
photometric redshifts (m« > 24) or even detections (m« > 25.5) 
from COMBO-17. Without a comprehensive multi-wavelength 
photometric redshift analysis we cannot definitively address 
this issue. Yet, we explore this issue briefly here to assess the 
potential scale of the biases, choosing to focus on two key is- 
sues: optical detectability and redshift accuracy. 

One key point is that of optical brightness. In terms of af- 
fecting the total IR luminosity density or skewing the fractions 
of IR luminosity from each morphological type, the most IR- 
luminous galaxies are the most relevant; i.e., luminous (10 12 Lq > 
L > 10 n L Q ; LIRGs) and ultra-luminous (> 1O 12 L ; ULIRGs) 
infrared galaxies. Thus, it is important to check if local ULIRGs 
and LIRGs would be visible if placed at z ~ 0.7. We compare 
the optical properties of a sample of local universe ULIRGs and 
LIRGs from Surace, Sanders, & Evans (2000) and Arribas et al. 
(2004) with the sample of Pec/lnt galaxies from the z ~ 0.7 slice 
(Fig. A7). Photometry for these galaxies was presented in B and 
/ bands; we corrected these for galactic foreground extinction, 
and transformed them to U and V absolute magnitudes using 
the PEGASE stellar population model (if instead one assumed 
a very young underlying stellar population, and postulated that 
redder colors were due to dust alone, the transformation from 
B—I to U — V would remain unchanged to < 0.1 mag, owing to 
the strength of the age/metallicity/dust degeneracy in the opti- 
cal wavelength region). In spite of their high obscurations, local 
ULIRGs and LIRGs are reasonably bright in the optical, and 
in common with z ~ 0.7 Pec/lnt galaxies, show a wide range 
in galaxy colors (primarily because of the dust obscuration). 
Highly-obscured IR-luminous galaxies are not faint in the op- 
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Az > 0.2 is 7/142, or 5% 22 . The RMS and outlier fraction do 
not strongly depend on IR luminosity. 

Thus, bearing in mind that z ~ 0.7 IR-luminous galaxies 
are reasonably bright in the optical, and that COMBO- 17's 
classifier appears to work equally well for IR-bright and IR- 
faint galaxies, we tentatively conclude that the philosophy used 
in this paper — exploring the IR properties of an optically- 
selected sample in a thin redshift slice — should not lead to any 
serious biases against IR-luminous, heavily obscured galaxies. 
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FIG. A8. — A comparison of VVDS spectroscopic and COMBO-17 pho- 
tometric redshifts as a function of 24^im flux. Diamonds denote high-quality 
spectroscopic redshifts and points denote those with more uncertain line iden- 
tifications. 199 galaxies with My < —20 and 0.6 < z < 0.85 are shown. The 
grey vertical line approximately shows the canonical 10 u Lq limit for LIRGs, 
showing a considerable LIRG population in this field and redshift range. Ar- 
rows denote 24fim upper limits. 



tical; they are rather fainter than they would be in the absence 
of dust, but they are still reasonably bright, and are certainly 
bright enough to be easily picked up at z ~ 0.7 to an ras-band 
limit of 24 21 . 

Yet, it is not clear that these dusty galaxies will be well- 
classified by COMBO-17's photometric redshift classifier. COMBO- 
17 adopts a 2-dimensional stellar population and SMC extinc- 
tion curve dust reddening grid for classification, allowing it in 
principle to classify galaxies with a wide diversity of stellar 
population properties and dust contents. In Fig. A8 we test the 
performance of COMBO-17's classifier directly using redshifts 
from the VIMOS VLT Deep Survey (VVDS; Le Fevre et al. 
2004) as a function of the 24/im flux. To ensure a large sam- 
ple, we expand the redshift range to 0.6 < z p hot < 0.85, and we 
require that My < -20, as this range contains the vast majority 
of all MIPS detections. There are a total of 199 galaxies with 
VVDS redshifts with estimated photometric redshifts and rest- 
frame V-band magnitudes in this range. Diamonds show galax- 
ies with well-measured spectroscopic redshifts (quality flags 
3/4/23/24), whereas lower-quality spectroscopic redshifts, typ- 
ically with less certain line identifications, are shown as points 
(quality flags 1/2/9/21/22). It is clear that most galaxies have 
well-constrained photometric redshifts; the robust sigma of the 
distribution is <7 m \, usl (Az) = 0.04, independent of IR luminos- 
ity. Furthermore, the fraction of LIRGs with Az > 0.2 is 5/57, 
or 9%, whereas the fraction of galaxies with L < 1O U L0 with 

21 ULIRGs and LIRGs at z > 1 will have m R > 24; therefore, optically- 
selected redshift surveys to this limit will be unable to construct unbiased sam- 
ples of IR-luminous galaxies. 



B. A COMPARISON OF GOODS AND GEMS MORPHOLOGIES 
OF IR-BRIGHT GALAXIES 

An unavoidable source of error comes from classification un- 
certainties. With the classifications based on the 1 -orbit rest- 
frame V-band data from GEMS, > 70% of the SFR is in undis- 
turbed galaxies. Yet, when one probes to fainter surface bright- 
ness levels, one will be able to better recognize tidal tails, lead- 
ing to an increased fraction of interacting galaxies. To esti- 
mate the importance of this source of uncertainty, we classified 
a subset of 290 0.65 < z < 0.75 galaxies with deeper F850LP 
data from GOODS (Giavalisco et al. 2004, the GOODS south- 
ern dataset forms the central <~ 1/4 of the GEMS area, giving 
F850LP data with 5x longer exposure time). Overall, around 
twice as many galaxies appeared to be clearly interacting in 
the deeper dataset; the detailed increase in fraction varied from 
classifier to classifier. Almost all of the newly-recognized in- 
teractions were relatively faint — GEMS -depth data are insuf- 
ficient to unambiguously discern the faint tidal features that are 
visible in the deeper GOODS data. This issue is discussed at 
length inWolfetal. (2005). 

Owing to small number statistics, it is impractical to sim- 
ply analyze the properties of the GOODS sample alone. We 
therefore derive a correction to the IR luminosity function that 
is tuned to reproduce the systematic difference between the 
GOODS and GEMS classifications for the 290-galaxy subsam- 
ple covered by GOODS. For each galaxy, the nearest 25 galax- 
ies in terms of IR luminosity are selected. Then relative split 
of spiral : E/S0 : Irr : Pec/Int for these 25 galaxies using 
the GEMS classifications is derived, and is compared with the 
type split derived using the GOODS classifications. This pro- 
cess is repeated for all galaxies, eventually building up the split 
of galaxy types as a function of IR luminosity in GEMS and 
GOODS. The difference between these type splits, as a function 
of IR luminosity, is used to correct the GEMS type-dependent 
IR luminosity functions in Fig. 3 to roughly compensate for the 
shallower depth of GEMS compared to GOODS. The main dif- 
ferences between GOODS and GEMS is that GOODS-depth 
data shows a statistically-significant increased fraction of clear 
interactions at log 10 L IR /L Q ~ 11, while the fraction of irregu- 
lar and spiral galaxies at these luminosities is slightly reduced 
to compensate. 

This issue will be discussed further by Papovich et al. (in 
preparation). 

22 The worst outliers have lower quality VVDS spectroscopic redshifts; Wolf 
(priv. comm.) finds in a more comprehensive comparison that the VVDS spec- 
troscopic redshifts with low confidence quality flags 1/2/9/21/22 have a ~ 25% 
error rate, and that COMBO-17 has an overall error rate of ~ 5% or less. 
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